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ABSTRACT 

Redshift surveys like the Sloan Digital Sky Survey (SDSS) have given a very precise 
measurement of the galaxy luminosity function down to about Mr = —17 (« Mb = 
— 16). Fainter absolute magnitudes cannot be probed because of the flux limit required 
for spectroscopy. Wide-field surveys of nearby groups using mosaic CCDs on large 
telescopes are able to reach much fainter absolute magnitudes, about Mr = —10. 
These diffuse, spiral-rich groups are thought to be typical environments for galaxies 
so their luminosity functions should be the same as the field luminosity function. 
The luminosity function of the groups at the bright end {Mr < —17) is limited by 
Poisson statistics and is far less precise than that derived from redshift surveys. Here 
we combine the results of the SDSS and the surveys of nearby groups and supplement 
the results with studies of Local Group galaxies in order to determine the galaxy 
luminosity function over the entire range —25 < Mr^ < —9. The average logarithmic 
slope of the field luminosity function between Mr = —19 and Mr = — 9 is a = —1.26, 
although a single power law is a poor fit to the data over the entire magnitude range. 
We also determine the luminosity function of galaxy clusters and demonstrate that 
it is different from the field luminosity function at a high level of significance: there 
are many more dwarf galaxies in clusters than in the field, due to a rise in the cluster 
luminosity function of a —1.6 between Mr = —17 and Mr — —14. 

Key words: galaxies: luminosity function, mass function 



1 INTRODUCTION 

The galaxy luminosty function (LF) (j}{L) is defined such 
that 4>{L) dL is the density of galaxies having luminosities 
between L and L -I- dL. It is usually regarded as a conve- 
nient test of galaxy formation theories because it is so easy 
to measure. However, the observed galaxy LF is difi'erent 
from the galaxy mass function, which is more directly pre- 
dicted by theory. This difi'erence is particularly significant 
in the context of cold dark matter galaxy formation theories 
(Moore et al. 1999; Klypin et al. 1999; Bullock, Kravstov & 
Weinberg 2001; Chiu, Gnedin & Ostriker 2001; Benson et 
al. 2003a,b; Kazantzidis et al. 2004). 

There are a number of redshift surveys that have re- 
cently been completed. These surveys measure redshifts and 
therefore distances for large numbers of galaxies in fiux- 
fimited samples. The largest such survey is the Sloan Digital 
Sky Survey (SDSS; www.sdss.org) which has measured red- 
shifts for > 10^ galaxies. The LF produced by this survey 
has very small uncertainties. The Poisson errors are negli- 
gible due to the sample size, and the sample is complete 
because the vast majority of high-luminosity galaxies have 



high surface brightnesses (Cross et al. 2001) and so are con- 
tained within the sample. A number of other surveys (see 
Section 2) also have produced accurate LFs, all of which are 
consistent with each other given the uncertainties that come 
from matching the different filters used in the different sur- 
veys and from cosmic variance. The SDSS redshift sample is 
only complete to absolute magnitudes of about Mr = —17.5. 
This limit is imposed by the requirement that galaxies need 
to be bright enough for spectroscopy. The spectroscopy limit 
is quite bright because the SDSS uses a relatively small (2.5 
m) telescope. 

Another recent development has been the advent 
of mosaic CCDs that can be used on large telescopes 
(e.g. SuprimeCam on the NAOJ Subaru 8 m Telescope; 
Miyazaki et al. 2002). This has permitted observations of 
large areas of sky down to very faint limits. A number 
of low-surface-brightness galaxies are seen in fields within 
galaxy groups but not in blank fields. These are interpreted 
as dwarf galaxies in the group, not background galaxies 
(Flint et al. 2001a). We can therefore infer their distances 
without spectroscopy, which is not possible for these faint 
low-surface-brightness galaxies. Hence measurements of the 
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LF in nearby (< 20 Mpc) groups are possible down to 
Mr = -10. 

Diffuse, spiral-rich groups are typical environments for 
most galaxies. Therefore it should be possible to append the 

LFs of those groups to the field LF from SDSS and obtain 
the galaxy LF over a very wide magnitude range. This is the 
subject of the present paper. 

At the very faint end {Mr ~ —10), the surveys of 
groups run into difficulty because of the inability to assign 
memberships to faint low-surfacc-briglitncss galaxies. At 
R > 22, low-surface-brightness {jj, > 25 R mag arcsec"^ over 
2 or more seeing disks) galaxies appear in blank field sur- 
veys. Thus, while excess numbers of low-surfacc-brightncss 
galaxies may be seen in the group fields, wc can only infer 
group membership in a statistical sense. For example, none 
of the candidate galaxies with Mr > —11 in the diffuse 
groups studied by Trentham & Tully (2002) can be regarded 
as highly probable members on an individual basis. 

The only environment where the LF can be measured 
fainter than Mr = —10 is the Local Group (the faintest 
galaxy in the Local Group, Andromeda IX, has Mr = —9; 
Zucker et al. 2004). We therefore supplement the combined 
SDSS-|-Groups LF with the Local Group LF at the very faint 
end. 

The paper is structured as follows. In Section 2 we 
describe results from redshift surveys. In Section 3 wc in- 
vestigate the behaviour of the LF at very bright magni- 
tudes {Mr < —23). In Section 4 we describe results from 
deep wide-field surveys of groups. In Section 5 we com- 
pute the i?-band Local Group LF. In Section 6 we combine 
all the individual LFs to determine the galaxy LF between 
Mr = —25 and Mr = —9. In Section 7 we investigate an- 
alytic forms that describe this function. In Section 8 wc do 
a brief compaxison with theory. In Section 9 we compute 
the cluster LF and investigate how it is different from the 
field LF. Finally in Section 10 wc summarize. Throughout 
this work, we assume the following cosmological parameters: 
Ho = 65 kms-^Mpc-\ ^matter = 0.3, Qa = 0.7. 



2 RESULTS FROM REDSHIFT SURVEYS 

The SDSS LF is most appropriate for this study because its 
r' filter is similar to the R filter that was used in CCD mosaic 
studies of nearby groups. It is derived from measurements 
of > W'' galaxies so the statistical uncertainties are small. 
There are a number of other redshift surveys which have 
similarly produced field LFs (Table 1). 

The SDSS LF is described in detail by Blanton et 
al. (2001, 2003). In Figure 1 we present the SDSS LF con- 
verted to the R magnitude system that we use in this paper. 
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Figure 1. The SDSS luminosity function, from Blanton et 
al. (2003), corrected to our cosmology. This is derived from the 
SDSS r' LF assuming a colour correction oi r' — R varying from 
0.24 at Mr = -24 to 0.17 at Mr = -16. These conversions 

come from the caliibrations of Fukugita, Shimasaku and Ichikawa 
(1995) and the colour-dependent LFs of Blanton et al. (2001). 
Three values of a = d log(ji(L)/d logL, the logarithmic slope of 
the luminosity function, are shown. 

There are too few galaxies in the local Universe {D < 40 
Mpc; Nearby Galaxies Catalog, Tully 1988) for this region of 
the LF to be measured directly: only one galaxy (the Seyfert 
galaxy NGC 4594) has Mr < -24 and 25 have Mr < -23. 
For D < 40 Mpc, the ratio of galaxies with -24 < Mr < 
-23.5 to those with -23.5 < Mr < -23 is 0.20 ± 0.11. 
The LF shown in Figure 1 implies a value for this ratio 
of 0.09, consistent with the number for the local Universe 
number. This concordance suggests that the evolutionary 
corrections made to the SDSS data, which are computed 
using observations of galaxies whose absolute magnitudes 
span a wide range, are appropriate at the very bright end 
Mr < -23. 



3 THE VERY BRIGHT END 

At the very bright end the SDSS LF is not a direct mea- 
surement. This is because most of the galaxies in the SDSS 
sample with Mr < —23 are sufficiently distant that evolu- 
tionary corrections to redshift zero are significant. Figure 2 
shows that at absolute magnitudes Mr < —23 the major- 
ity of galaxies have redshifts z > 0.2, where corrections are 
typically 0.3 mag (Lin et al. 1999; Blanton et al. 2003). 



4 RESULTS FROM DEEP MOSAIC CCD 
SURVEYS 

In Table 2 we list some nearby groups where deep CCD 
mosaic imaging has been performed. The LFs from these 
surveys extend down to Mr = —10. Earlier photographic 
studies (e. g. Ferguson & Sandage 1991) covered larger areas 
but did not reach magnitudes this faint. All of these groups 
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Table 1. Spectroscopic field surveys 



Survey 


Reference 


Filter 


Limiting Magnitude 


Stromolo-APM 


Loveday et al. 1992 


bj 


-17 


Hawaii-Caltech 


Cowie et al. 1996 


B 


-12.5 






K 


-20.5 


Autofib 


Ellis et al. 1996 


bj 


-15 


LCRS 


Lin et al. 1996 


R 


-18 


ESO Slice Project 


Zucca et al. 1997 


bj 


-15 


SDSS 


Blanton et al. 2001, 2003 


u' 


-16.5 






9' 


-17.5 






r' 


-17 






i' 


-18 






z> 


-19.5 


2MASS +2df 


Cole et al. 2001 


J 


-19 






Ks 


-20 


2df 


Norberg et al. 2002 


bj 


-17.5 




Figure 2. The fraction of galaxies at redshift z > 0.2 in the SDSS 
sample used to compile Figure 1. 



have boon imaged in the Cousins R band (Acft ~ 6588 A) , 
which is close to the SDSS r' band (Aeff = 6290 A; Fukugita 
et al. 1995). The 7?-band LF of these groups is therefore 
suitable for extending the SDSS LF to faint magnitudes. 

Members are identified based on their surface bright- 
nesses. Low-surfacc-brightncss galaxies that appear within 
a group of giant galaxies and whose counterparts do not ex- 
ist in blank fields are assumed to lie at the distance of the 
group. Distances are therefore determined based on mor- 
phology; spectroscopic redshifts are unattainable for these 
faint low-surfacc-brightncss galaxies. At the faintest magni- 
tudes this exercise becomes statistical as a few very small 
low-surface-brightness galaxies are visible in blank fields. We 
tend to see more such galaxies in fields containing nearby 
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Figure 3. The composite luminosity function for the groups 
listed in Table 2, with errors computed as decribed in the text. 
Here A^g^i is the total number of galaxies in the combined sample, 
proportional to N in Figure 1. 



groups of large galaxies, but for any particular galaxy we 
can only assign membership on probabilistic grounds. 

The composite LF of the groups in Table 2 is presented 
in Figure 3. Computing this was not straightforward be- 
cause of the different methods used by different authors to 
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Table 2. Groups with deep luminosity function measurements 



Group 


Distance 
Mpc 


limiting Mn 


reference 


survey telescope 


sample size 


Ursa Major 


18.6 


-11 


Trentham ct al. 2001a (TTVOl) 


CFHT 


50 


Leo I 


11.1 


-10 


Flint et al. 2001b 


KPNO 0.9 m 


112 


Coma I 


16.4 


-10 


Trentham & TuUy 2002 (TT02) 


Subaru 


38 


NGC 1023 Group 


10.0 


-10 


Trentham & Tully 2002 (TT02) 


Subaru 


28 



assess membership probabihties. For Ursa Major, Coma I, 
and the NGC 1023 Group (hereafter Gl), TTVOl and TT02 
employed a subjective rating system based on surface bright- 
nesses, morphologies and spectroscopy (when available). For 
the present calculation, galaxies were assigned a high (rat- 
ing "0-2" in the notation of TTVOl and TT02) or low (rat- 
ing "3") probability of membership. For the Leo I Group 
(hereafter G2), Flint, Bohe & Mendes de Olivera (2001b) 
computed completeness corrections C(M)in each magnitude 
interval. The number of galaxies in the composite group 
sample in each magnitude interval equals the number of Gl 
galaxies with a high membership probability plus 1/2 times 
the the number of Gl galaxies with a low membership prob- 
ability plus (1 -|- C(M)) times the number of galaxies in the 
G2 sample. The weighting of 1/2 for the Gl galaxies with 
a low membership probability comes from a study of the 
radial distribution of galaxies with this rating within the 
groups (TTVOl, TT02). 

The uncertainties in these numbers are potentially large 
because the rating system used by TTVOl and TT02 is sub- 
jective and the completeness corrections used for the G2 
sample are model-dependant. Errors in the number of group 
members will be systematic and a conservatve approach will 
be required to determine the uncertainties; we estimate them 
from the quadrature sum of Poisson errors (Gehrels 1986), 
the entire contribution of Gl galaxies with a low member- 
ship probability and the entire enhancement in the number 
of G2 galaxies due to completeness corrections. That the er- 
ror bars in Figure 3 are so large at the faintest magnitudes 
follows from the fact that the contributions from the last 
two terms dominate the uncertainties. 

The LF of the groups over the range Mr = —19 to 
Af_R — —10 is fit by a power law of index a = —1.26 ± 0.11. 
This is an interesting magnitude range because Mr = — 19 
{Mb ~ —18) is the brightest magnitude where dwarf galax- 
ies, as identified by their position on a magnitude-surface- 
brightness plot (e. g. Figure 1 of Binggeli 1994), contribute 
significantly to the LF. 

Precise distances have been measured to individual low- 
luminosity galaxies in Centaurus and Sculptor (.Jerjen, Free- 
man & Binggeli 2000; Karachentsev et al. 2002). These 
groups will be good candidates to add to the list in Ta- 
ble 2 when complete galaxy samples are available; this will 
reduce the errors in Figure 3. 



5 THE LOCAL GROUP 

The LF of groups has large errors at Mr ~ — 10 which can 
be decreased by including the Local Group in the analysis. 
In the Local Group, the lowest-luminosity galaxies can be 
resolved into stars so that membership can be established 
with confidence. However, there are very few galaxies and 
the Poisson errors are large (but not significantly larger than 
the errors in the LF of groups at A^r ~ —10). 

The faintest galaxies in the Local Group would have 
been detected in the nearby groups studied by Trentham & 
Tully 2002 given the detection limit of that survey {la of 
about 29 R mag arcsee^^). However only a small part of 
the galaxies would be visible above the sky and the galax- 
ies would be indistinguishable from a large number of small 
background galaxies, many of which have been dimmed by 
(l-l-z)* cosmological effects. The ability to distinguish mem- 
bers from background galaxies can only be performed with 
any confidence brighter than Mr = —10. 

There are, however, complications in deriving a Local 
Group LF that can used in the current analysis. 

Firstly there are uncertainties as to the membership of 
stars in specific galaxies. For example, in Draco, inclusion 
of a number of new stars caused the tidal radius to increase 
significantly (Piatek et al. 2001 Odenkircheu et al. 2001) 
which in turn caused the absolute magnitude to brighten. 
Therefore absolute magnitudes of individual Local Group 
galaxies are uncertain. 

Secondly, most Local Group galaxies have not been ob- 
served in the i?-band. Therefore we need to perform a colour 
correction for each dwarf. 

Thirdly, the Local Group sample may be incomplete at 
the faint end. The lowest- luminosity galaxies have very low 
surface brightnesses and so can be difficult to find. They are 
particularly difficult to find if they are viewed in projection 
against either the Milky Way or the M31 disks (e. g. An- 
dromeda IX; Zucker et al. 2004). They are also difficult to 
find if they are in the process of being disrupted (e. g. Sagit- 
tarius; Ibata, Gilmore & Irwin 1994, 1995). 

Furthermore, it can be possible to identifty groups of 
stars but be uncertain as to whether or not these stars form 
part of a discrete galaxy with a dark halo. Including or 
excluding such objects from the Local Group galaxy sam- 
ple can then be subjective. An example is Andromeda VIII 
(Morrison et al. 2003), which may be part of the M31 stream 
(Ibata et al. 2004). Another example is Canis Major (Martin 
et al. 2004, Momany et al. 2004). 

In Table 3 we present a list of Local Group members and 
their absolute Ji-band magnitudes. The methods used to de- 




Figure 4. The luminosity function of the Local Group. 



rive these magnitudes arc: [1] interpolation using the spec- 
tral energy distributions (SEDs) of Fukugita ct al. (1995; 
hereafter F95), adopting a literature value for the absolute 
magnitude in some filter, and choosing a SED appropriate 
for the Hubble type of the galaxy; [2] as [1] but choosing a 
SED appropriate for a galaxy with broadband colours equal 
to those measured for the galaxy; [3] as [1] but choosing a 
SED appropriate for a gala^xy with broadband colours equal 
to those inferred from colour-magnitude diagrams; [4] adopt- 
ing a literature value for Mr. 

The LF for the Local Group is presented in Figure 4. 
The faint-ond slope isa= —1.1 ±0.1, consistent with the 
findings of Pritchet & van den Bergh (1999) and van den 
Bergh (2000). 

Recently Karachcntscv ct al. (2004) compiled a list of 
galaxies in the local Universe with distances less than 10 
Mpc, the Local Volume (LV) catalogue. This sample is prob- 
ably incomplete at the very faint end at present but over the 
next few years it should be possible to identify and measure 
magnitudes and distances (by, for example, the tip of the 
red giant branch method) for the vast majority of galaxies 
within 10 Mpc. This will allow us to generate a LF with 
much smaller errors than in Figure 4. 



6 COMBINING THE MEASUREMENTS - THE 
GALAXY LF 
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Table 3. Local Group members 



Galaxy 


type 


Mb 


method 


rererence 


M31 = NGC 224 


Sb 


-21.8 


1 


1 


Milky Way 


Sbc 


—21.5 


1 


2 


M33 = NGC 598 


So 


— 19.5 


3 


2,3 


LMC 


dl 


— 18.8 


3 


2 4 


SMC 


dl 


— 17.4 


3 


2 5 


M32 = NGC 221 


E 


—17.1 


2 


6 7 


NGC 205 


dE 


—16.9 


2 


6 7 


NGC 6822 


dl 


—16.5 


2 


6 7 


IC 10 


dl 


— 16.3 


2 


2 8 


NGC 3109 


dl 


— 16.1 


2 


6,7 


NGC 185 


dE 


-15.7 


2 


6,7 


IC 1613 


dl 


— 15.8 


3 


6 9 


NGC 147 


dE 


— 15.7 


2 


6,7 


Sagittarius 


dE 


—15.4 


3 


6,10 


Sextans A 


dl 


—14.9 


3 


6,11 


WLM 


dl 


—14.7 


3 


6,12 




dl 


—14.7 


3 


6,13 


Forii&x 


dE 


— 13.6 


3 


6 14 




dl 


— 13.4 


3 


6 15 


And VII — Cassiopeia 


dE 


— 12.7 


3 


6,16 


Leo I 


dE 


— 12.6 


3 


6 17 


And I 


dE 


— 12.5 


3 


6 18 


And II 


dE 


—12.5 


3 


6,18 


SaffDIG 


dE 


—12.3 


4 


19 


Leo A 


dl 


—12.1 


3 


6,20 


Antlia 


dE 


—11.9 


3 


6,21 


Aquarius = DDO 210 


dE 


— 11.5 


3 


6 22 


And VI = Pegasus 


dE 


-11.1 


3 


6,23 


LGS 3 


dl 


-10.9 


3 


24,25 


And III 


dE 


-10.8 


3 


18,26 


Cctus 


dE 


-10.7 


3 


6,27 


Leo II 


dE 


-10.6 


3 


6,28 


Sculptor 


dE 


-10.3 


3 


6,29 


Phoenix 


dE/I 


-10.3 


3 


6,30 


Tucana 


dE 


-10.2 


3 


6,31 


Sextans 


dE 


-10.1 


3 


26,32 


Draco 


dE 


-10.0 


3 


6,33 


And V 


dE 


-9.8 


3 


34,35 


Carina 


dE 


-9.9 


3 


1,11,15 


Ursa Minor 


dE 


-9.8 


3 


34,35 


And IX 


dE 


-8.8 


3 


1,10 



1. SED from NASA/IPAC Extragalactic Database (ned- 
www.ipac.caltech.edu); 2. van den Bergh 2000; 3. McConnachie 
ct al. 2004; 4. Castro ct al. 2001; 5. Alciano ct al. 2003; 6. Grcbcl 
ct al. 2003; 7. Prugnicl & Hcraudcau 1998; 8. de Vaucoulcurs 
& Abies 1965; 9. Tikhonov & Galazutdinova 2002; 10. Ibata 
et al. 1997; 11. Dolphin et al. 2003; 12. Rejkuba et al. 2000; 
13. Sakai et al. 1997; 14. Bersier & Wood 2002; 15. Gallagher 
et al. 1998; 16. Grebel & Guhatakurta 1999; 17. Held et al. 2001; 
18. Da Costa, Armandroflf & Caldwell 2002; 19. Lee & Kim 2000; 
20. Dolphin ct al. 2002; 21. Tolstoy & Irwin 2000; 22. Lee et 
al. 1999; 23. Hopp et al. 1999; 24. Lee 1995a; 25. Miller et al. 2001; 
26. Caldwell ct al. 1992; 27. Sarajadini ct al. 2001; 28. Lcc 1995b; 
29. Monkicwicz ct al. 1999; 30. Gallart ct al. 2001; 31. CastcUani, 
Marconi & Buonanno 2002; 32. BcUazzini ct al. 2001; 33. Klcsscn, 
Grebel & Harbeck 2003; 34. Caldwell 1999; 35. Davidge et 
al. 2002; 36. Rizzi et al. 2003; 37. Zucker et al. 2004 



The field luminosity function, obtained by combining the 
LFs described in the previous three sections, is presented in 
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Figure 5. The field galaxy luminosity function. The normal- 
ization is as appropriate to the SDSS data and the other sam- 
ples are grafted onto the SDSS LF. The dashed line shows 
the best Schechter function fit to all the data: = —22.0, 

a* = —1.29. The dottcd-dashed line shows the best power-law fit 
for Mji > -19; a = -1.26. 

Figure 5 and Table 4. Prior to being combined, the LFs were 
normalized to a consistent scale and were weighted at each 
magnitude by the inverse of the square of the uncertainty. 

Figure 6 shows how much each individual LF con- 
tributes to the total LF. The SDSS dominates the LF at 
the bright end and the CCD mosaic surveys at the faint 
end. The Local Group is important only for the faintest 
{Mr = -9.5) point. 

Two points on the LF are worthy of further discussion: 
(1) the Mr = —17.5 point. This is the faintest point 
where the SDSS dominates and it is significantly lower than 
all other points in this region. One possibility is that the 
SDSS sample is incomplete here beacuse galaxies with this 
absolute magnitude have surface brightnesses too low for 
spectroscopy. On the other hand, the majority of galax- 
ies with Mr = —17.5 {Mb ~ —16.5) have surface bright- 
nesses within their half-light-radius of /iso ~ 20 i? mag 
arcsec^^ (Trentham, TuUy & Verheijen 2001b), which is 
bright enough for spectroscopy with the instruments used 
by SDSS. Another possibility is that galaxies are rare at 
precisely this magnitude, which correponds roughly to the 
transition magnitude between giant and dwarf galaxies. In- 
terestngly, Flint, Bolte & Mendes de Olivera (2003) report a 
deficiency in galaxies with absolute magnitudes close to this 
value in the Leo I Group. 



Table 4. The Field Galaxy Luminosity Function 



Mr 4> 

N mag~i Mpc~* 



-25.5 (2.18 ± 2.18) x 10"*^ 

-25.0 (2.86 ± 0.46) X 10-9 

-24.5 (3.17 ± 1.64) x 10"* 

-24.0 (7.26 ± 1.40) x lO"'' 

-23.5 (1.30 ± 0.12) X 10-5 

-23.0 (8.38 ± 0.36) x lO'^ 

-22.5 (3.50 ± 0.091) x 10"* 

-22.0 0.000920 ± 0.000026 

-21.5 0.00164 ± 0.000052 

-21.0 0.00224 ± 0.000092 

-20.5 0.00288 ± 0.00015 

-20.0 0.00285 ± 0.000096 

-19.5 0.00339 ± 0.00014 

-19.0 0.00390 ± 0.00019 

-18.5 0.00483 ± 0.00035 

-18.0 0.00573 ± 0.00072 

-17.5 0.00325 ± 0.00075 

-16.5 0.00694 ± 0.0019 

-15.5 0.00931 ± 0.0022 

-14.5 0.0150 ±0.0029 

-13.5 0.0156 ±0.0028 

-12.5 0.0228 ± 0.0037 

-11.5 0.0166 ± 0.0033 

-10.5 0.0393 ± 0.0083 

-9.5 0.0151 ± 0.0073 



T — • — ' — ' — • — I — • — ' — ' — • — I — • — ' — ' — • — r 




-25 -20 -15 -10 

Figure 6. The weighting of each component to the total LF 
for the SDSS (solid line), CCD mosaic (dashed line) and Local 
Group (dotted-dashed line) samples. The weightings here are pro- 
portional to (T^, where a is the error in the component LF. 
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(2) the Mr — —9.5 point. There are only three Local Group 
galcixies which contribute to this point, and all narrowly 
escaped inclusion into the Mr = —10.5 bin. While com- 
pleteness and small number statistics are surely concerns at 
this extreme faint end, this observation may be a hint of 
a turnover in the LF at about Mr = —9.5. Galaxies like 
Andromeda IX would then be extremely rare. The shape of 
the LF here and the possible existence of a turnover should 
become well-established over the next few years as the LV 
sample is extended to larger distances and completeness is- 
sues in that sample become better understood. 

The normalization of the LF is affected by cosmic vari- 
ance and will vary from survey to survey. The value adopted 
in Figure 5 comes from SDSS observations which cover 0.3 % 
of the sky and so should provide a reasonable representation 
of the cosmological average. 



7 ANALYTIC FORMS 

A power-law fit to the LF over the range Mr > — 19 gives a 
logarithmic slope a = —1.2 but this is not a good fit in that 
it overpredicts the number of galaxies at the extreme faint 
end. The error bars in Figure 5 are small enough that cur- 
vature in the LF within this magnitude range is detectable 
at a high level of significance. 

It is common to fit LFs with Schechter (1976) functions: 

0(L) = 0*exp(-^)(^)" i_, (1) 

or in magnitude units 

0(M) = 0.92 r (iol-o-<^-^*)])-*+i ^-lol-o-t"--*)! _ 

(2) 

Here (f>* is a normalization density, L* and M* = 
—2.5 logj^Q L* arc a characteristic luminosity and magnitude, 
and a* is a characteristic faint-end slope. These fits normally 
have a > —2 and so provide convenient ways to describe 
galaxy luminosities over a magnitude range (close to M*) 
where most of the total luminosity of a sample of galaxies 
is obtained. Schechter function fits have no physical basis 
but are attractive in that their exponential cutoff is simi- 
lar to that predicted by simple analytical models of galaxy 
formation (e. g. White & Rees 1978). 

A Schechter function fit over the entire range —25 < 
Mr < —9 is poor because (1) it is not possible to simul- 
taneously fit the data at both Mr = —22 and Mr = —20 
where the errors in the SDSS LF are so small, (2) the point at 
Mr = —17.5 is systematically too low, and (3) a Schechter 
function which fits the data well at bright magnitudes then 
overpredicts the number of faint galax;ies. These features are 
illustrated by the dashed line on Figure 5. 

Two-parameter fits cannot fit the data over any signif- 
icant magnitude range. Composite forms, like linear combi- 
nations of Gaussian and Schechter functions, may provide 
acceptable fits to the data but there is no physical basis for 
selecting a particular set of functions at this stage. 



8 COMPARISON WITH THEORY 

It is well known that cold dark matter (CDM) theory pre- 
dicts a mass function of dark halos that is steep at the low- 



mass end (e. g. Moore et al. 1999). Simulations predict a 
logarithmic mass function slope a ~ — 2 at low masses. 

The field galaxy LF is much shallower than this at all 
magnitudes. This tells us that star formation must be less 
efficient in lower-mass galaxies. There are various physical 
processes that can be responsible for this. Lower-mass gala:x- 
ies are more susceptible to gas ejection by supernovardriven 
winds (Dekel & Silk 1986) due to their having smaller grav- 
itational potential wells. This means that negative feedback 
effects during star formation (Efstathiou 2000) arc stronger 
in lower-mass galaxies. Additionally, in hierarchical cluster- 
ing cosmologies like CDM, many low-mass galaxies form at 
later times so it will bo difficult for them to collect gas from 
an intergalactic medium that is being progressively heated 
by the metagalactic ultraviolet background (Tlioul & Wein- 
berg 1996; Klypin et al. 1999; Bullock, Kravstov & Weinberg 
2001; Tully et al. 2002). 

Whatever processes are responsible, we shall see below 
that they must operate less efficiently in galaxy clusters. 
Perhaps the presence of a confining medium (Babul & Rees 
1992) inhibits the effects of negative feedback during low- 
mass galaxy formation (see Section 5 of Roberts et al. 2004). 



9 THE CLUSTER LF 

Various measurements of cluster LFs have been made pos- 
sible by a number of advances in technology. 

The prototypical rich galax;y cluster is the Coma cluster, 
with several hundred galax;ies brighter than Mr = — 19 per 
square degree. Recently a spectroscopic LF of Coma extend- 
ing down to Mr = — 16 has been published (Mobasher et 
al. 2003). This is a difficult measurement to make because 
it requires spectroscopy of large numbers of faint sources, 
many of which have relatively low surface brightnesses. This 
measurement has been made possible by the advent of wide- 
field multi-object spectrographs that can be used on large 
telescopes, in this case WYFFOS on the 4.2 m William Her- 
schel Telescope (Bridges 1998). 

Many other clusters that are richer and more distant 
than Coma have also been studied (A665 at z = 0.18: the 
only AbeU (1958) Richness 5 cluster, Wilson et al. 1997, 
Trentham 1998; A963 at ^; = 0.21, Driver et al. 1994, Tren- 
tham 1998; A1689 at z = 0.18, Wilson et al. 1997; A868 
at z = 0.15, Boyce et al. 2001; A1146 at z = 0.14, Tren- 
tham 1997). These clusters have very well-determined LFs 
at the bright end because their contrast against the back- 
ground field galaxy population is so high. Photometry alone 
can therefore be used for a measurement of the LF at bright 
absolute magnitudes. When many such clusters are consid- 
ered in conjunction with each other, the numbers of galax:ies 
in each bin becomes large and the Poisson errors are conse- 
quently small. At the faint end, these measurements cannot 
be used to determine the LF since the cluster counts are 
small compared to the background counts and the field-to- 
field variance of the background is much larger than the 
cluster signal. Spectroscopic measurements cannot help ei- 
ther at the faint end since the clusters are so distant. 

The Virgo cluster is much less rich than the Coma Clus- 
ter, but at 17 Mpc it is sufficiently nearby that dwarf galaxies 
are much larger than the seeing and can be identified by their 
low surface brightnesses. The LF can then be determined 
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to faint limits by photometric observations alone. This LF 
will not include the contribution from compact galaxies like 
M32 and the Ultra Compact Dwarfs in the Fornax Clus- 
ter (Phillipps et al. 2001, Drinkwater et al. 2003) because 
they will be mis-idcntificd as background ellipticals. How- 
ever these types of galaxies are rare. The Virgo Cluster Cat- 
alog (Binggeli, Sandage & Tammann 1985) was compiled 
about 20 years ago and provided a precise measurement of 
the bright end of Virgo Cluster LF (Sandage, Binggeli & 
Tammann 1985). More recently, mosaic CCD measurements 
(Trentham & Hodgkin 2002, Sabatini et al. 2003) using the 
Wide Field Camera on the 2.5 m Isaac Newton Telescope 
(INT) have confirmed the earlier results and have extended 
the LF fairitward by about 5 magnitudes. Extremely deep 
observations over a smaller area of Virgo using SuprimeCam 
on the 8 m Subaru Telescope (Trentham & Tully 2002) did 
not result in the discovery of any new dwarfs brighter than 
the completeness limit of the INT survey (about Mb = — 11) 
that were too low in surface brightnesses to be seen with the 
smaller telescope. The implication is then that the LF of the 
Virgo Cluster is known at least as faint as Mb = — 11. 

We derive a composite cluster LF by combining mea- 
surements from photometric studies of rich galaxy clusters 
(A665, A963, A1795, A2199 and A1146; Trentham 1997, 
1998), the deep spectroscopic study of the Coma cluster 
(Mobasher et al. 2003), and a photometric study of the Virgo 
cluster (Trentham & Hodgkin 2002). Each dataset is impor- 
tant in a different absolute magnitude regime (Figure 7). 
At very bright absolute magnitudes, the rich cluster sam- 
ple is most important; the Coma and Virgo samples suffer 
from small number statistics. At intermediate magnitudes, 
the Coma sample is most important; the rich cluster LF 
has large errors since a background subtraction is required 
and the background has a large field-to-field variance and 
the Virgo sample has substantial Poisson errors. At faint 
magnitudes only the Virgo data can be used; the other sur- 
veys cannot detect dwarfs, most of which have low surface 
brightnesses. Although the three samples utilise very differ- 
ent methods in establishing membership, they are consistent 
with each other in the absolute magnitude where they over- 
lap. All data used are _R-band data except the Virgo data 
which are B-band data converted to the _R-band using a 
combination of (1) B — R colours from a calibration de- 
rived using measurements of faint galaxies observed in both 
the INT and Subaru surveys, (2) the spectral energy distri- 
butions from F95, and (3) the B — R calibration for giant 
galaxies of different Hubble Type T obtained using data 
from Tully & Pierce (2000): B - R = lAQ - 0.059T (equa- 
tion (1) in TT02). Typical random errors in this conversion 
are estimated to be 0.1 - 0.2 mag, smaller than the bin size 
of 0.5 mag. 

The cluster LF is presented in Table 5 and Figure 8. 
This LF is appropriate for centres of galaxy clusters al- 
though the inclusion or exclusion of cD galaxies in the sam- 
ples used to define the LF is arbitrary. There is evidence, 
however, that the LF might be quite different in the outer 
parts of galaxy clusters (Phillipps et al. 1998). 

The cluster LF can be fit by a Schechter function bright- 
ward of Mr = —16, where it rises with a logarithmic slope 
a = —1.6 and then is flat faintward of about Mr = —13. 
The faintest point on the LF shown in Figure 8 might be af- 
fected by incompleteness (although the deep Subaru Virgo 
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Figure 7. The weighting of each component to the total LF. 
The lines are for the rich cluster (solid line). Coma (dashed line) 
and Virgo (dotted-dashed line) samples. The weightings here are 
proportional to iT~^, where a is the error in the component LF. 



Table 5. The Galaxy Cluster Luminosity Function 



Mr 4> 

N mag~i Mpc~^ 



-23.0 


0.88 it 0.61 


-22.5 


1.46 ±0.80 


-22.0 


16.2 ±4.0 


-21.5 


33.4 ±5.8 


-21.0 


57.3 ±7.2 


-20.5 


78.5 ± 8.6 


-20.0 


84.2 ±9.2 


-19.5 


72.9 ±8.9 


-19.0 


89.1 ± 10.4 


-18.5 


113 ± 12 


-18.0 


101 ± 13 


-17.5 


152 ± 11 


-17.0 


150 ± 14 


-16.5 


214 ± 19 


-16.0 


180 ± 26 


-15.5 


366 ± 65 


-15.0 


404 ± 71 


-14.5 


599 ± 89 


-14.0 


551 ± 87 


-13.5 


577 ± 89 


-13.0 


561 ± 88 


-12.5 


558 ± 88 


-12.0 


460 ± 80 


-11.5 


250 ± 59 
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Figure 9. The logarithmic slope a of the field (open circles) 
and cluster (filled circles) LFs, and the difi^erence between the 
two (open triangles). The points are calculated at each absolute 
magnitude Mjj by considering the shape of the LF over a two 
magnitude interval centred on M^. 



Figure 8. The galaxy cluster luminosity function, computed as 
described as in the text. The normalization is arbritrarily chosen 
to be that appropriate for the Coma Cluster at Mji = —21. The 
dashed line is the best-fitting field luminosity function, which we 



approximate by a Sclicchtcr function with Mt 



-22.0 and q* = 



-1.28 brightward of Mr 
faintwajrd of Mjj = —19. 



-19 and a power law with a = —1.24. 



measurements suggest otherwise) but the flattening faint- 
ward of about Mr = — 13 comes from an analysis of the 
Virgo data in an absolute magnitude regime where the com- 
pleteness is close to 100 per cent. 

The field and cluster LFs are highly inconsistent with 
each other (Figure 9). Brighter than Mr = -20, the LF falls 
off slightly more steeply in the field than in clusters due to 
the presence of luminous elliptical galaxies in clusters. At 
about Mr = —20, the field and cluster LFs have similar 
shapes, despite being dominated by contributions from dif- 
ferent kinds of galaxies (early-type galaxies in clusters and 
late-type galaxies in the field) . The cluster LF becomes sig- 
nificantly steeper than the field LF at about Mr = — 17 due 
to the presence of a large number of dwarf elliptical galaxies. 
Both LFs are reasonably flat faintward of about Mr = —15, 
although the cluster LF has a higher normalization relative 
to the giant galaxies than does the field LF and the contri- 
bution of dwarf irregular galaxies relative to that of dwarf 
elliptical galaxies is higher in the field than in clusters. This 
excess number of cluster relative to field dwarfs is described 
and studied in detail by Roberts ct al. (2004), although the 
differential between the cluster and field that we find is less 
than that found by those authors. This may be due in part 
due to there being slightly fewer dwarfs per giant in groups 



than in environments with very low galaxy densities (see 
Zabludoff & Mulchaey 2000). 

The faint end of the cluster luminosity function depends 
on measurements of the Virgo Cluster and it is important 
to establish that this environment is representative - differ- 
ent studies of the faint end of the Virgo LF (Trentham & 
Hodgkin 2002, Trentham & Tully 2002, Sabatini ct al. 2003) 
concur as regards the general shape (rise at Mr = —17 
and flattening at fainter magnitudes) but comparisons be- 
tweeen different clusters at the very faint end have yet to 
be made. There is a hint of a LF with a similar shape in 
the knot of early-type galaxies around NGC 1407 (Tren- 
tham & Tully 2002) but this sample is small. The best place 
to determine a LF that can be compared with the Virgo 
LF is the Fornax Cluster, which also has a population of 
low-surface-brightness dwarf galaxies (Kambas et al. 2000). 
Spectroscopic information is available for a subset (about 
one-quarter) of the Fornax; Cluster Spectroscopic Survey 
(FCSS) dwarfs, but there are too few dwarfs known (19 with 
-15.8 < Mb < -12.7; Deady et al. 2002) for a LF to be 
computed that can be compared in detail with the Virgo 
LF. When the FCSS is complete, this will be an important 
comparison. 

The previous paragraph is concerned with the precise 
shape of the LF at the faintest magnitudes. All current ev- 
idence does at this stage points towards an excess of low- 
luminosity galaxies in clusters, although this precise shape 
depends on the measurements of the local clusters. For ex- 
ample, an excess of low luminosity cluster galaxies above 
the expected number given the field LF was seen in the 2dF 
redshift survey (de Propris et al. 2003). 

This tendency for dwarfs to be more numerous in the 
richer and denser environments is important in that it tells 
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us that the physical processes responsible for suppressing 
the formation of stars in low-mass dwarf galaxies may op- 
erate less efficiently in clusters than in the field. Alterna- 
tively there are additional physical processees that operate 
preferentially in clusters and alleviate the effects of dwarf 
suppression there. For a more detailed discussion of these 
issues, the reader is referred to Roberts et al. (2004). 

The shape of the cluster LF is such that it cannot 
be described by a simple analytic form. Composite forms 
may be more appropriate; these may provide a useful test 
of theories which may wish to invoke different formation 
mechanisms for giant and dwarf galaxies. For example, 
the cluster LF is well fit by a double Schcchter function 
with (A/*, a*) = (-21.4,-1.0) for giants and (M*,q*) = 
(—16.6,-1.1) for dwarfs and a relative normalization of 

2.7 X 10-°-*(*^iant-*^d„arf). 



10 SUMMARY 

This paper presents the results of a study that combines 
different datasets in order to measure the field galaxy LF 
over the range —25 < AIr < —9. The data come from the 
Sloan Digital Sky Survey, from CCD mosaic imaging surveys 
of nearby groups, and from the Local Group. The field lu- 
minosity function is well-described by a Schechter function 
with Mr = -22.0 and a* = -1.28 brightward of Mr ^ -19 
and by a power law with a — —1.24 faintward of Mr = —19. 
This is shallower than the CDM mass function (a ~ —2), 
implying that star formation is less efficient in lower-mass 
galaxies. 

A similar analysis was performed to measure the LF of 
galaxy clusters. The three samples used were a photometric 
survey of rich Abell clusters, a spectroscopic survey of the 
Coma Cluster and a photometric survey of the Virgo Clus- 
ter. The cluster LF differs from the field LF in that it has 
a rise of a ~ —1.6 between Mr = —17 and Mr = —14, 
although it fiattens faintward of Mr = —14. In the context 
of CDM theories, this suggests that the physical processes 
that are responsible for suppressing the formation of stars 
in low-mass galaxies operate less efficiently in clusters than 
in the field. 
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